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We present in detail a nonperturbative method for vortex liquid systems. This method is based 
on the resummation of an infinite subset of Feynman diagrams, the so-called parquet graphs, con- 
tributing to the four-point vertex function of the Ginzburg-Landau model for a superconductor in a 
magnetic field. We derive a set of coupled integral equations, the parquet equations, governing the 
structure factor of the two-dimensional vortex liquid system with and without random impurities 
and the three-dimensional system in the absence of disorder. For the pure two-dimensional system, 
we simplify the parquet equations considerably and obtain one simple equation for the structure 
factor. In two dimensions, we solve the parquet equations numerically and find growing translational 
order characterized by a length scale R c as the temperature is lowered. The temperature depen- 
dence of R c is obtained in both pure and weakly disordered cases. The effect of disorder appears as 
a smooth decrease of R c as the strength of disorder increases. 

PACS: 74.20.De,74.60.Ge 



I. INTRODUCTION 

Since the discovery of high- T c superconductors, the na- 
ture of the mixed state in a type II superconductor has 
been a focus of theoretical and experimental investiga- 
tion. In many high-T c materials, thermal fluctuations 
are responsible for the melting of the vortex lattice phase 
predicted by the mean field theoryB into a vortex liquid 
phase. The presence of quenched random impurities in 
the vortex liquid or lattice phases also plays an important 
role, since it presents a possibility of a dissipation-free 
current flow due to pinning of flux lines. New theoretical 
phases such as the vortexn-glass phascB for point defects 
and the Bose glass phaseH in the presence of extended 
defects have been-,proposed. As shown by Larkin and 
Ovchinnikov (LO)q the quenched point disorder destroys, 
for spatial dimension d < 4, the long-range crystalline 
order of the vortex lattice. The system is described by 
some characteristic length scale R c over which a short- 
range translational order exists. 

In a previous papeiQ, we developed a nonperturba- 
tive scheme to calculate the structure factor of the two- 
dimensional (2D) vortex liquid in the absence of random 
impurities. The main ingredient of this nonperturba- 
tive method was the resummation of an infinite subset 
of Feynman diagrams contributing to the structure fac- 
tor, the so-called parquet graphsB. This is an analytic 
approach to 2D vortex liquid system which is sophisti- 
cated enough to predict growing crystalline order in the 
system as the temperature is lowered. The growth of the 
translational order was investigated in connection with 
the sharp peaks developing in the liquid structure fac- 
tor. Within this scheme, we found no evidence for a fi- 
nite temperature phase transition into the vortex lattice 
phase and the system remains as a liquid. The length 
scale, R c , characterizing this growing translational order 
seemed to diverge only in the zero temperature limit. 

In this paper, we give detailed derivations of the par- 



quet graph resummation technique which were omitted 
in the previous papertj. In addition, we present a new 
simplified version of the parquet equations governing the 
structure factor, which were given in terms of a set of 
coupled nonlinear integral equations previously. In the 
present work, we were able to obtain one simple equa- 
tion for the structure factor, which contains all the non- 
perturbative information for the 2D vortex liquid in the 
absence of disorder. 

In this paper, we also apply the parquet graph resum- 
mation technique to the 2D vortex liquid in the presence 
of quenched disorder. We find that the sharp peaks which 
appeared in the structure factor of the pure system be- 
come broadened as the strength of the disorder increases. 
This is entirely consistent with one's intuition that in the 
presence of disorder the length scale R c describing the 
translational order becomes smaller as compared to the 
pure case. In the present work, since our pure system is 
always in the liquid phase, the effect of disorder appears 
as smooth deviations from the pure case. From the non- 
perturbative results, we find the temperature dependence 
of the length scale R c . It is, however, difficult to make 
any connection between our results and the LO type ar- 
gument, since there exists no vortex lattice phase at any 
finite temperature within our nonperturbative scheme, 
while the LO argument always starts from the vortex 
lattice with a true long-range crystalline order. 

We note that there exist recent theoretical studiesQ 
based on the elastic theory of pinned lattices suggest- 
ing that the pinning by quenched disorder become less 
effective due to the periodicity of the lattice so that a 
quasi long-range order persist beyond the Larkin length 
scale R c . In the present nonperturbative analysis, the ex- 
istence of such quasi long-range translational order has 
not been observed. 

The point whether the 2D vortex liquid in the absence 
of disorder undergoes a finite temperature phase transi- 
tion into a 2D vortex lattice is still controversial. Nu- 
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merical simulations! seem to suggest a first-order phase 
transition. However, as shown in Ref. |j, these results 
depend cruicially on boundary conditions. In a spherical 
geometry, the authors of Ref. [9] demonstrated the absence 
of a finite temperature ukase transition. There also ex- 
ists a recent experimentE3 performed on a sample with 
very weak pinning, where no sign of a phase transition 
is detected. The present parquet approximation, which 
is an analytic theory on an infinite plane, seems to sup- 
port the absence of a finite temperature phase transition 
between 2D vortex liquid and solid. 

Unlike the 2D system, it is generally believed that a 
3D vortex liquid undergoes a first-order phase transition 
into a vortex lattice. An experiment by Zeldov et al.t3 
on BSCCO has been accepted as a convincing evidence 
for the transition. However, there exists a recent claiint 2 ] 
that the results of Ref. [ll] might be an artifact due to a 
particular sample geometry. Therefore, it would be inter- 
esting if one could develop a three-dimensional parquet 
graph resummation scheme and obtain nonperturbative 
information on the 3D vortex liquid system. We find 
that one can generalize the parquet equations to three 
dimensions without difficulty. Unfortunately, the equa- 
tions become very complicated and we were not able to 
obtain a numerical solution to the 3D parquet equations. 

In the next section, we introduce the structure factor 
of a two-dimensional disordered vortex liquid within the 
Ginzburg-Landau theory. In Sees. [II and [Tv|, we present 
detailed derivations of the parquet equations which ac- 
count for all the parquet graphs contributing to the struc- 
ture factor for both pure and disordered cases. We also 
consider zero-dimensional models to discuss the validity 
of the parquet approximation in general. In the following 
section, we present the main results of our calculation and 
discuss the temperature dependence of the length scale 
R c . In Sec. VI, we present the generalization of the par- 



quet equations to three dimensions. Finally, we conclude 
with discussion on future directions. 



II. THE STRUCTURE FACTOR 

We begin our analysis with the Ginzburg-Landau free 
energy for a superconducting film in a perpendicular 
magnetic field B = V x A in the presence of quenched 
random impurities; 

M*l = [ d 2 r( — \(-mV - e*A)4-| 2 
J \2m 

+ (a + T (r))|v|>| 2 + ||*| 4 ), (2.1) 

where a, (3 and m are phenomenological parameters. The 
random field r(r) representing the quenched impurities 
satisfies the probability distribution, r(r) = and 

r(r)r(r') = A<5 (2) (r - r'). (2.2) 



In this paper we neglect the fluctuations in the vector po- 
tential A and restrict the order parameter 'J to the space 
spanned by the lowest Landau level (LLL) wavefunctions. 
In the symmetric gauge, where A = y,x), the LLL 
is fully described by an arbitrary analytic function of the 
variable z — x + iy multiplied by an exponential factor; 
*(x,y) = exp(-^z*z)0(z), where \j? = e*B/fi = 2w/Q 
and Q is the area of the unit cell of the vortex lattice. In 
the LLL approximation, the GL free energy becomes 



F [(f)] = J dz*dzf(a H +T(z,z*)) exp(-^-| 



+ ^exp(V N 2 M ^)| 4 



(2.3) 



where J dz*dz denotes the integration over x-y plane and 
an = a + he* B /2m vanishes at the mean field transition 
temperature. 

The central quantity in this analysis is the structure 
factor of the two-dimensional vortex liquid. It is pro- 
portional to the Fourier transform of the density-density 
correlation function, x(k) = J d 2 R e lk R x(r, r + R); 



X (r,r') = (|*(r)H*(r')| 2 ) - <|*(r)P) (|*(r')| 2 ), 

where the angular brackets denote the thermal averages. 
The structure factor A(k) is then defined by 



Qexp( 



2/i ; 



■)A(k)=x(k)/ 



-I 2 



|*(r)|* 



(2.4) 



A convenient way to deal with quenched averages is to 
introduce n replicas of Z and calculate the correlation 
functions with respect to 



Z n = Yl ex P 



dz*dz 



{ Qff ^I^V|4(z)| 2 



"M \ z \ 



i^ Q (z)n<Mzr} 



(2.5) 



in the limit n — > 0. On e can develop a standard per- 
turbation theory for (2.5). The bare propagator arising 
from the perturbation expansion of (2.5) is given byt-3: 



= o ab — — exp(— ( z), 
an Zir 2 



(2.6) 



where the double bracket ((. . .)) denotes the average with 
respect to Z n . It is important to note that, because of 
the simple quadratic term of the LLL free energy in ( [2.5| ), 
the renormalized propagator is obtained by simply 
replacing an by the renormalized mass an. This means 
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that the magnetic length is the only length scale in- 
volved in the renormalized propagator. This is one of 
the simplifying features of the LLL approximation that 
makes the present nonperturbative calculation feasible. 
As noticed in Ref. O, the GL free energy obtained from 



(2.5) is not closed under renormalization. The renormal- 
ization effectively generates the quartic vertices of a gen- 
eral form; 



j~J dz*dzie ■ 



Ml^il'+M 2 ) 



(2.7) 



i=l,2 



53 |0a(^i)| 2 (s ab g(\ Zl - z 2 \) - w(\zi - z 2 \)\ \Mz2)\ 2 ] , 

- a,b ^ ' * 



for arbitrary functions g and w. The bare interactions 
correspond to 

9b(\z\) = ^5(z)5(z*), w b (\z\) = ~5(z)5(z*). 

It is convenient to work with the Fourier transformQ 

g(k) = J dz*dz g(\z\)exp(^(k*z + kz*j^ 

and similarly u>(k), where fc, k* are complex momenta, 
k = ki + ik 2 ,k* = ki — ik 2 with k = {ki,k 2 ). Thus, 
9s{\k\) = [3/2 and ws(|fc|) = A/2 are constants. 



In order to calculate the structure factor, (2^), we 
need to consider the renor mal ized four-point correlation 
function arising from Eq. (2.5); 



m(4W b (4)M^)M^))) = G a R c (zl,z 3 )G b *(z* 2 ,z 4 ) 
+G%(z*,z 3 )G a R d (zl,z 4 ) + ({PM4>tte)M*a)M**)))c, 

where the last term denotes the contribution from all 
connected Feynman diagrams. To the lowest order of 
perturbation theory, the connected four-point correlation 
function can easily be calculated as 

m(zm(4)M^)M^)))c = -<*a C <W^r(g)) 2 



a 



A 1 ( z l Z3+Z2 2<l)/2 



xe 



dk*dk 



5 a bgB(\k\) - w B (\k\) 



(27T) 2 

i(k'(z 3 — z 4 )+k(zl — z^))/2 

~d, z 3 ^ z 4 ) + 0(/3 2 ,/?A,A 2 ), 



(2.8) 



where the second term on the right hand side is the same 
as the first t erm with c and d, and z 3 and 2:4 inter- 
changed. In (2.8), we absorbed the renormalization of 



the propagators using the renormalized mass a R . To 
the low est order, one can easily evaluate the k-integrals 
in ( |2.8[ ), but, in general, higher order correction s to the 
connected four-point function are represented in (2.8) by 
the departure of the quartic vertex functions from con- 
stants, gs(|fc|) = (3/2 and WB{\k\) = A/2, to general 
k-dependent functions, g R (k) and w R (k). Therefore, 



the LLL approximation enables us to concentrate on the 
renormalized quartic vertex function g R (k) and w R (k), 
which depend only on one variable instead of three char- 
acterizing three independent channels in a usual field the- 
ory. 



Now, one can put (2^) for general <?i?(k) and w R (k) 
in a more symmetric form using scaled functions, 

/ JJ (k) = |exp(-k 2 /2 A1 2 ).9 i? (k), 

P 

v R (k) = ^exp(-k 2 /2/i 2 )^(k). 

The bare vertices are given by /s(fc) = vs{k) = 
exp(— k 2 /2^t 2 ). Interchanging z% and Z4 in the second 
term on the right hand side of (2^) is equivalent to using, 
instead of f R (k) and v R (k), the transformed functions, 
/^(k) and w}f{(k)_where /(k) is defined for an arbitrary 
function /(k) byEl 



/(k) 



2tt 
7- 



d 2 _p 
d 2 p 



1 

(2tt) 2 ^( p ) exp (^2^ lP2 ~ fc2Pl ))' 



7(p) exp ( — (kiP2 - k 2 pi)). 



The general four-point function is then given by (see Ap- 
pendix) 



■ 2 - <■ dk*dk /3 



IA2 



(2.9) 



x ex P ( y^ Z3 + Z * Z4 ' 



(2tt 



2 2 -^afc,cd(k) 



exp ( - ^2 - ~(k*(z 3 - 24) + fc(z* - 4)j), 



where 



r a 6,cd(k) = -S ac S bd S ab f R (k) - 9v R (k) 



<Wfl(k) - 0v R (k) 



(2.10) 



and 9 = A//3 represents the strength of the disorder. 
Eq. (|2. 10 ) can be represented diagrammatically as in 
Fig. \f~^ 

The structure factor A(k) is then obtained by joining 
t wo external legs of the four -point correlation fu nctio ns, 
(pL9l ). From and ( |2.10| ) and the definition Q) in 

the limit n —>■ 0, we obtain 



where 

r(k) 

and x 
which ap; 



A(k) = 1 - 2xr(k), (2.11) 



f R (k) + f R (k) - 8 [v R (k) + v R (k)] (2.12) 



- " 2 /?/2 7 ra 2 f 

s in the high-temeprature perturbation 



[i 2 f3/2ira R is a dimensionless parameter 



expansion^ 



3 




a, 3 a,1 a, 3 




a, 2 a, 4 b,2 




l-x(2- 6)+x- 



,2tt 



d 2 k 



2(1 - 0)/ fl (k) - 0(2 - 0)v fl (k) 



(2-13) 



where olt = a//^/27r//3^ 2 is the dimensionless temper- 
ature. Note that this is an exact relationship between 
the renormalized propagator and the quartic vertex func- 
tions. We will, however, use /i?(k) and fi?(k) obtained 
from the present nonperturbative approximation. The 
Hartree approximation used in the high-temperature per- 
turbation expansion!!^ corresponds to neglecting terms 
that depend on the vertex functions, fu and vr. 



FIG. 1. Diagrammatic representation of the renormal- 
ized connected four-point correlation function. The labels, 
1, •••,4, denote z*,z%,zs,Z4, respectively, and a, ■•■,d are 
replica indices. 
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FIG. 2. Diagrammatic representation of Dyson equation 
for the self-energy E. The solid lines represent the renor- 
malized propagators. The thick dashed and dotted lines are 
the renormalized vertex functions, fa and vr, as in Fig. 1, 
while the thin dashed and dotted lines denote the bare vertex 
functions, Jb and vb, respectively. 

In the following sections, we will evaluate /i?(k) 
and vji(h) nonperturbatively by summing the parquet 
graphs. As mentioned earlier, we absorb any renor- 
malization of the propagator into the renormalized pa- 
rameter aR. But knowledge of the four-point vertex 
function also fixes the relationship between an and the 
bare parameter thus completing the description of 
the system. This relation comes from the Dyson equa- 
tion arising from ( |2.5[ ) which is described diagrammat- 
ically in Fig. @ for the self-energy X = G _1 — Gq 1 — 



III. PARQUET GRAPH RESUMMATION: PURE 
SYSTEM (A = 0) 

In order to calculate //s(k) and f/s(k), one needs to 
evaluate the Feynman diagrams contributing to the four- 
point correlation function. In this and the next sections, 
we show in detail how to sum an infinite subset of such 
diagrams, the so-called parquet graphs. 

We first consider the pure system where A = 0, then 
generalize the result to the disordered case. The analysis 
of the pure case will mostly reproduce results given in 
Ref. |^. In the present work, we make a further simplifi- 
cation of the parquet equation for /i?(k) found in Ref. |B| 
and obtain a very simple equation for the structure factor 
A(k). 

For both pure and disordered cases, we will present 
the analysis for the zero-dimensional analogues of (^), 
which can be integrated exactly. For the pure case, we 
find it convenient to introduce the parquet resummation 
scheme in a simple zero-dimensional model and then gen- 
eralize to the two-dimensional problem. Furthermore, 
since there is no apparent expansion parameter involved 
in the parquet approximation, it would be instructive to 
apply the parquet resummation to the cases where ex- 
act solutions are known and to compare the approximate 
result with the exact solution. 



A. d=0 

For A = and d ime nsion d = 0, the partition function 
corresponding to (2.3) is a simple integral, 



dipdip* 2 /3. , 4 . 

— exp{-a H \ip\ - -\ip\ ) 

Z7T Z 

2^ erfc(cr) exp(cr 2 ), 



where a = an/ 'V^P and erfc(er) = 1 — crf(cr) is the com- 
plementary error function. In the parquet analysis, one 
calculates the renormalized four-point vertex 
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where cv^ 1 = (|V'| 2 )- Although an exact solution can 
readily be found, one can construct the usual Feynman 
graph expansion for T. To the lowest order, we have 
r = Tb = (3/2. To the one-loop order, the diagrams can 
be represented as shown in Fig. |^. A convenient way to 
generate the next order diagrams is to replace each ver- 
tex in a one-loop diagram by the vertices obtained up to 
the one-loop order. In general, one can construct higher- 
order diagrams by replacing each vertex in the one-loop 
diagrams by the vertices obtained up to the current order 
of perturbation expansion. An example of such construc- 
tion is shown in Fig. |3|. The diagrams obtained in this 
way are called parquet graphs. Note that parquet graphs 
can be separated into two parts by cutting two propaga- 
tor lines. 



;xx: :>x; :x< 



(a) 



>X>><x< 
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+ o(/f) 
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FIG. 4. The totally irreducible vertex R. 





FIG. 5. The parquet decomposition of reducible vertex IT. 
Again, the labels 1, • • • 4 are drawn for the general cases. The 
same diagrammatic decomposition can be used in higher di- 
mensional models for both pure and disordered cases. 



(b) 

FIG. 3. (a) Three one-loop diagrams. The symmetry 
factors are 2,4 and 4, respectively. The labels 1, • • • 4 are 
drawn for the general cases to be discussed later. For the 
zero-dimensional case, there is no distinction between the sec- 
ond and the third diagrams. But, it is important to separate 
this contribution into two equal parts, (b) An example of 
constructing successive parquet graphs. 

Although parquet graphs cover an enormous number of 
diagrams, obviously not all diagrams can be constructed 
in this way. The non-parquet diagrams are generated in 
the above construction by the so-called totally irreducible 
vertices whose contribution here is denoted by R. The 
totally irreducible vertex consists of the bare vertex and 
higher order (0(/3 4 )) vertices (see Fig. ^). There is no 
systematic ways of enumerating these higher order di- 
agrams contributing to R. The parquet approximation 
which we employ in this work corresponds to keeping only 
the bare vertex contribution in R: R ~ Tb = (3/2 in the 
zero-dimensional case. 



If both bare vertices in a one-loop diagram are replaced 
by the full renormalized vertex in the above construction 
of graphs, then the diagrams become overcounted. There 
is a systematic way to eliminate this overcounting and to 
generate all the Feynman diagrams and the symmetry 
factors associated with each diagram, and that is to use 
the full vertex for the bare vertex on the right hand side 
of Fig. ^|, but the so-called irreducible vertex for the one 
on the left hand side, which is defined as follows: Let 
IT, i — 1,2,3 denote the diagrams constructed out of the 
three one-loop diagrams in Fig. || and Aj the correspond- 
ing irreducible vertices, we have from Fig. |5| 

2 

Hi = T A 1 T, (3.1a) 

a R 

n 2 = -4rA 2 r, (3.1b) 

a R 

n 3 = -^A 3 r, (3.ic) 

a R 

where we absorb any renormalization on propagator lines 
into aR. Now, for the irreducible vertex Ai, we must in- 
clude all the diagrams that do not belong to IT to avoid 



the overcounting, thus, Aj = T — Ilj. Finally, the renor- 
malized vertex T is given by the sum of all contributions; 



Therefore, 



r = r 



A, = R 



5> 



En. 



(3.2) 



(3.3) 



Note that (3.1)-(3.3) are exact relations for the renormal- 
ized vertex T. But, as mentioned before, we use R ~ j3/2. 

One can easily sim plify the above equations for R = 
(3/2 as follows: From (3.1), one has 



III 



2r 2 



2r' 



n, = n. 



AT 2 



a|-4r- 



Therefore, from (3.2), we obtain 

„2 27 2 



7 



7 



1-7 1-27 



(3.4) 



where 7 = 2T /a 2 R and p = (3/a 2 R . As a function of p, 
there is only one solution j(p) that satisfies the trivial 
condition: T = 7 = when (3 = p = 0. 




FIG. 6. The renormalized four-point vertex function 7 for 
the pure zero-dimensional model as a function of tempera- 
ture a. The solid line is the exact solution. The dashed line 
corresponds to the parquet approximation. 

To complete the description of the system, one needs 
a relation between the bare (an or a) and renormalized 
(an or p) mass. In the parquet approximation, we use an 
analogue of Dyson equation, ( 2.13j ), which, in this case, 
is 



1 



l-2p(l- 7 ) 



(3.5) 



From this relation using 7 = ^(p) obtained from (|3.4|), 
one finds p as a function of <r, and consequently we have 
the renormalized four-point vertex 7 as a function of a. 
In Fig. H 7(c) of the parquet approximation is compared 
with the exact solution. One can see there is excellent 
agreement between the two. 



B. d=2 

We now consider the two-dimensional problem in the 
absence of disorder. The parquet equations can be con - 
stru cted similarly to the d = case. In fact, eqs. (3.1)- 
( |3.3|) take almost the same form as before. The only dif- 
ferences are the fact that the vertices now depend o n the 
internal momentum k and the right hand sides of (3.1) 
should be obtained from a direct evaluation of diagrams 
in Fig. |^. In Appendix, we show explicitly how this cal- 
culation is done. The resulting equations corresponding 
to (IO) are 



Hi(k) 

n 2 (k) 
n 3 (k) 



-x(Aior)(k), 

-2xA 2 (k)r(k), 
-2x(A 3 *r)(k), 



(3.6a) 
(3.6b) 
(3.6c) 



where the operations o and * between two arbitrary func- 
tions /(k) and g(k) are defined by 



(/o fl )(k) = H£ 



(/*0)(k) 



2tt 
,2 



d 2 p 

d 2 p 



/(k - p)s(p) cos ( 



/(k-p)g(p). 



hP2 - k 2 pi 



A* 



FJ (2tt) 2 

The remaining parquet equations are given by 

3 

r(k) = i?(k) + En i (k), 

i=l 

A i (k) = J R(k) + En J -(k), 

3& 



(3.7) 
(3.8) 



where i?(k) represents the totally irreducible part. In the 
parquet approximation, R(k) is equal to the bare vertex 
part: i?(k) ~ (l/2)(/ B + f B ) = cx P (-k 2 /2p 2 ). 

We can make a simplification on the parquet equations 
as in the zero-dimensional case. Wejirst note the follow- 
ing identities: (fog) = (fog), fog = f o g = fog, 

and / * 5 (k) = /(k)g(k). Inserting A, (k) = r(k)-IT(k) 
into ( 3.6b| ) and ( 3.6c ) and using r = T from ( 2.1 0| ) , we 
obtain 



n 2 (k) = n 3 (k) 



-2a;r 2 (k) 
1 - 2xT(k) ' 



(3.9) 



Now, inserting this into (3.6a) 
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ni(k) = -xUr + n 2 + n 3 ) o rj (k) 

= -x((R + 2Tl 2 ) or)(k). (3.10) 



Therefore, inserting (3.9) and (3.1C) into (3.7), we obtain 
an equation for T(k). In terms of the structure factor 
A(k) = 1 — 2xr(k), this equation becomes 



1 - A(k) 
A(k) 



xR{k) +x{Ro A)(k) 
■(1-A) 2 



(3.11) 



A (k). 



This is our main equation mentioned in Sec. | that com- 
pletely describes the structure factor of 2D vortex liq- 
uid in the absence of disorder. Note that we have kept 
the totally irreducible part 7?(k) in a general form. We 
note that, although this is an exact relations in a very 
simple form^-it has little advantages for finding numeri- 
cal solutionsta over the coupled parquet equations, (3.6), 
(3.7) and (3.8|) or another version to be described be- 



low. We believe, however, that this equation might open 
a possibility in the future for a nonperturbative analytic 
investigation on the 2D vortex liquid. 

When we consider the disordered case in the next sec- 
tion, it will be necessary to decompose T(k) into /ii(k) 
and fn(k) using ( 2.10D and consider the parquet equation 
in terms of /j?(k) arising from the faithf ul r epresentation 
of Feynman diagrams as in Fig. [I]. Eq. (3.9) suggests the 
following decompositions: 

n 1 (k) = i(r 1 (k) + f 1 (k)), 
n 2 (k)=n 3 (k) = i(r 2 (k) + r 3 (k)), 

A 1 (k) = i(/ 1 (k)+7 1 (k)), 

A 2 (k) = A 3 (k) = i(7 2 (k) + 7 3 (k)), 



for some functions 1% and Ti. Inserting these into (3.6) 
we have 



r 1 (k) = -x(7 1 o/ fl )(k) (3.12a) 
T 2 (k) = -x(l 2 (k)f R (k)+I 2 (k)f R (k) 



and 



+ 7 3 (k)/ fl (k) 
r 8 (k) = -a:(/3*/fl)(k), 



/ B (k) = / fl (k) + ^r < (k) J 



(3.12b) 
(3.12c) 

(3.12d) 



Ji(k) = / B (k)+]Tr,(k), (3.12e) 

This version of the parquet equations has been previously 
given in Ref. [| 



IV. PARQUET GRAPH RESUMMATION: 
DISORDERED CASE 



For the disordered case, one has to construct the 
parquet equation for the vertex function containing 
n— replica indices. Let us first consider the two- 
dimensional case directly. It is straightforward to gen- 
eralize (3.6) to the present case. By putting the replica 
indices in the diagrams in Fig. ^, we have 



nS od (k) = -^(Ag e/ or e/>od )(k), 
ni'Uk) = -2^A« c/ (k)r /6;e ,(k), 
ng od ( k ) = - 2x E ( A 2/ d * r Ace) (k). 



(4.1a) 
(4.1b) 
(4.1c) 



The remaining equations follow from ( p.7[ ) and (3 

3 

TabMV = Rab.M + E n aLrf( k )' 



Aitd( k ) = ^( k )+E n Sc d (k), 



(4.1d) 
(4.1e) 



where 



R ab , cd (k) ~ r^ cd (k) = 5 ac S bd (S ab - 6) exp(-k 2 /2 A1 2 ) 

in the parquet approximation. 

As one can see from (2.10), in order to take n — > 
limit, one needs to decompose the above equations as in 
t he p revious section and to get equations analogous to 
( |3.12 ). First, we note that 



Aitc(k) - A^ cd (k), n^ >dc (k) = n^ cd (k). 

Therefore, we can write, for some functions Ti, Sj, 7j and 

Ji, 



,(2 



(3) 



t(2) 



n£ od (k) = -6 ac S bd £ ai A(k) - 05! (k) 



1 



ng cd (k) = n^(k) = ^S ac S bd S ab T 2 (k) - 9E 2 (k) 



K 3 ) 



S ab T 3 {k) - 6Z 3 {k) 



and 



ASc rf ( k ) = ^ ac 5 bd [5abh(k) - 9,h{k) 



Aitrf(k) = AaMc(k) - §*«c*m 5 ab I 2 {k) - 9J 2 {k) 



+ 2$ad5 bc 

TO) 



S ab h(k) - 6Ji(k) 
1 



+-^s ad s bc 



5 ab I 3 (k)-6J 3 (k) 
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Inserting these expressions and (2.10) into (4.1) and 
eliminating one term that contains an explicit factor of 
n, we get the following set of parquet equations in the 
presence of disorder: 

r x (k) = -x{h o f R - 6h o v R - 8Ji o / a )(k), 

r 2 (k) = -x(/ 2 (k)/ fl (k) + / 2 (k) [/ fl (k) - ^n(k)' 

+ [j 3 (k)-0J 3 (k)]/ Jl (k)), 
r 3 (k) = -x(/ 3 * - 0/ 3 * wji - 0J 3 * /fl)(k), 



and 



Si(k) = a;fl(Jiow fl )(k), 

H 2 (k) = -x(/ 2 (k) Ui? (k) + J 2 (k)/ R (k) 

+J 2 (k) [/ fl (k) - 0u fl (k)" 



+ 



J 3 (k)-0J 3 (k) « K (k) 



S 3 (k) = ^(J 3 *^)(k) 



with 



/ R (k) = / B (k) + ^r i (k), 

3 

Vfl(k) =w B (k) + ^S,(k), 



and 



j i (k) = / B (k)+X)r i (k), 



J i (k)= WB (k) + ^H j (k). 



(4.2a) 
(4.2b) 

(4.3a) 
(4.3b) 



Before solving these complicated equations, we con- 
sider first the zero-dimensional toy model where the par- 
qu et eq uations reduce to a set of algebraic equations as 
in (O). 



A. d=0 

For d = 0, the partition function is again a simple 
integral, 

Z(a„,r,l3) = J ^f- exp(-(o„+r)|^f-|w 4 ) 



Averaging over r with respect to (2.2), we have 



-lnZ((r,/9,0) = -ln/3-a 2 -- 



(4.4) 



da'e- (CT - ,T ' )2/9 lnerfc(a'), 



where tr = a^/y^, and 9 = X/ j3 as before. 

The quantities corresponding to /«(k) and Uij(k) in 
the zero-dimensional case are defined by 



./it 2/? 



/3 



d_ 

dj3 



-2<| 

a 2 



|2\2 



InZ, 



v4 r 



/3 (3 



d(3 



where a^ 1 = (|i/>| 2 ) = — <91nZ/<9a#. These quantities 
can easily be evaluated as functions of er and using 

On the other hand, the parquet eq ua tion s for / B and 
vr take exactly the same form as (4.2), (4.3), if we make 
the following simplifications: the binary operations o 
and * become multiplications between two factors; the 
transformation / has no effect, / = /; x is replaced by 
p = @/ctR\ and finally / B = Vb = 1. It is then straight- 
forward to eliminate Ii (or 1^) and Ji (or S,-) from the 
parquet equations in favor of f' R = pfu and v' R = pvR. 
The parquet equations then reduce to the following cou- 
pled algebraic equations. 



Vr 



R 



(1 



f R + v' R )(l-2f' R + v' R ) 



Vr 



R 



2v' 



i + v' R (i-2f R + v' R y 



2v' R = dp. 



(4.5) 



(4.6) 



For given 9 and p, one can find the corresponding so- 
lutions f' R and v' R t o thes e equations. We can then use 
the Dyson equation (2.13); 



1 



1 - p(2 - 9) + P [2(f R - v' R ) - 9(2f' R - v' R )} 



to find a relation between p and a: p = p(a) for given 9. 
Using this, we obtain f' R (a) and v R (a) for given 9. These 
are shown in Fig. ^ together with the exact solutions. For 
9 < 1, the parquet results show excellent agreement with 
the exact solutions. As 9 becomes larger, however, the 
discrepancies between two solutions grow. This analysis 
for the zero-dimensional model suggests that the parquet 
approximation is in general very good when the strength 
of disorder is moderate. But, the diagrams omitted in 
the parquet approximation might produce quantitative 
errors in the strong disorder regime. We note, how- 
ever, that physical quantities in this approximation re- 
main smooth-as functions of 9 unlike other approxima- 
tion methodala on this system. 
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FIG. 7. The renormalized four-point vertex functions, f' R 
and v' R for the disordered zero-dimensional model as functions 
of temperature a. The filled triangles, circles and squares are 
obtained from the parquet approximation for 6 = 0.1, 1.0 and 
5.0, respectively. The solid lines are the corresponding exact 
solutions. 



V. RESULTS AND ANALYSIS 

For the case of the two-dimensional system, one has 
to solve a s et o f coup led integral equations for /ij(k) 
and Uij(k), (4.2) and (4.3), containing two parameters 
x and 9. We consider a rotationally symmetric case 
where all vertex functions depend only on K = |k|//x. 
We apply a similar numerical technique to the one used 
m Rcf. |. In practice, we find it convenient to work with 
h R {K) = f R {K) - 9v R (K) and v R (K). For given set 
of irreducible parts, {Ii(K)} and {Ji(K)}, eqs. (4J2) are 
coupled linear integral equations for h R and v R . We first 
solve these equations for h R for fixed v R by numerically 
inverting a matrixrj, and then solve for v R using the 
solution h R . We then update the irreducible parts using 
(4.3). The solution to the parquet equations are obtained 
by iterating this procedure. 

A fast convergence can be obtained if we choose the 
initial functions, {ii}, {Ji}, and v R close to the actual 
solutions. At high enough temperatures, it is sufficient to 
start from Ii — Ji — v R — exp(— K 2 /2). As one goes into 
the low temperature regime, it is necessary to use the so- 
lution at a temperature close to the desired temperature 
as initial functions. We face the same numerical difficulty 
as in Ref. || as the temperature is lowered, namely one 
has to use a finer mesh in K space as well as a larger 
cutoff in order to get a low temperature solution. In this 
case, we have a coupled set of equations, which requires 
an additional computing time. The minimum temper- 
ature we used was ar — —8.5 where the cutoff was at 
K = 15 and number of mesh points was 600. We also 
find it difficult to solve the parquet equations directly for 



arbitrarily large values of 9. Again, one needs to start 
from the actual solution for 9 close to the value for which 
the structure factor is to be calculated. In the present 
analysis, we were able to obtain A (A') for three values of 
9; 9 = 0.1,0.2 and 0.3. 

In Fig. @, we first present the renormalized propagator 
yfx ~ as a function of temperature «r for three val- 
ues of 9. Compared to the pure (9 = 0) case, ar(x) as 
a function of x shows very little deviations. In general, 
one finds that the same value of parameter x represents 
a slightly higher temperature as 9 increases. 



-9 



-7 



-6 

<y.j 



-3 



FIG. 8. The renormalized propagator (fj, 2 /2nan ~ y/x) as 
a function of temperature ax in d — 2. The filled squares are 
obtained from the parquet approximation in the pure case. 
The open triangles, filled triangles and open squares corre- 
spond to the parquet approximation for 6 = 0.1,0.2 and 0.3, 
respectively. 

In Figs. H and [To], the structure factor is plotted for var- 
ious values of ar and 9. The structure factor develops a 
collection of peaks around the reciprocal lattice vectors 
(RLV) of the triangular lattice. As the temperature is 
lowered for fixed 9, one can clearly see from Fig. |^ that 
the first peak grows with decreasing width. This can be 
interpreted as a growing short-range translational order 
in the disordered vortex liquid. The length scale R c over 
which this order exists can be obtained from the inverse 
width |-nf the first peak, or equivalently from the peak 
hcighttm. For fixed temperatures, we find that the length 
scale decreases as the strength of disorder increases. (See 
Fig. |l^.) The result is consistent with an intuitive pic- 
ture where the disorder prevents ordering on long length 
scales. In Fig. |ll|, we show the height of the first peak as a 
function of temperature for 9 = 0.1, 0.2 and 0.3. We find 
that the length scale R c grows as \ax\ with decreasing 
proportionality constant for increasing 9. We recall that, 
in the absence of disorder, the length scale characterizing 
a growing crystalline order also grows as \ar\- But, from 
Fig. [ll], we find that the rate of this growth gets smaller 







as the disorder gets stronger. In the next subsection, we 
show that one can derive this behavior analytically from 
the parquet equations. 



Usin g thes e fun ctions and the original parquet equations, 
(4.2) and (13), one can find a set of equations for T(k) 
as follows: 



A 




K 



FIG. 9. The structure factor for fixed 9 = 0.2 and for var- 
ious temperatures (dashed line: ar = —3.78, dotted line: 
ax = —5.34, dot-dashed line: olt = —6.51 and solid line: 
olt = —8.36). The arrows indicate the positions of RLV of 
the triangular lattice. Since the second and the third RLV 
are closely spaced, our solution could not resolve these peaks. 



A. Parquet Equations for 9 1 

For very weak disorder, 9 -C 1, one can extract from 
the parquet equations the small-# behavior of the struc- 
ture factor A(k) around the pure (9 = 0) solution. 
One can in fact derive the temperature dependence of 
the length scale R c for small 9 from the parquet equa- 
tions. In order to do that, we nee d t o find the par- 
quet equation analogous to fl3.6| ), (3.7) and (3.8) for 
r(k) = (l/2)Uk) + h R (k)}. Out of the functions used 
in (|4.2|) and (4.3), we define 



Ai(k) 
A 2 (k) 



7 2 (k)+7 3 (k)- 



(Ji(k) + Jx(k) 
J 2 (k) + J 3 (k) 



A 3 (k) = A 2 (k). 



and 



ni(k) 
n 2 (k) 



ri(k) + ri(k)- 
r 2 (k) + f 3 (k)- 



(Si(k) + H!(k) 
~ 2 (k)+S 3 (k) 



Hi(k) 

n 2 (k) 



r(k) = r B (k) + ^n l (k), 
A 4 (k) = r(k)-n 4 (k), 

= -as(Aor)(k), 
= -2x 



A 2 (k)r(k) - -j 2 (k) VR (k) 



n 3 (k) = n 2 (k), 

where 



r B (k) = (1 - 9) exp(-k 2 /2/x 2 



(5.1a) 
(5.1b) 

(5.1c) 
(5.1d) 
(5.1e) 

(5.2) 



and we have used the above-mentioned identities involv- 
ing the binary operations, o and *. Note that these equa- 
tions a re a lmost t he s ame as those given in t he pu re case, 
( |3.6| ), (3.7) and (3JS). One extra term in ( 5. Id ) comes 
from the subtraction of the diagram containg an explicit 
factor of n. There is of course the second equation for 
vr, which couples nontrivially to these equations. 













T Tt 
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FIG. 10. The structure factor for fixed x = 30, (or slightly 
varying aj; —6.67 < «t < —6.45) and for various values of 
9: 9 = (thick line), 0.1,0.2 and 0.3. The peaks are getting 
smaller as the strength of disorder (9) increases. 

Note that the structure factor A(k) is just given by 
A = 1 — 2xT, and we can simplify the above equations as 
done in the pure ca se to get one equation for the s truc- 
ture factor. From (5. Id) and J 2 equation in (4.3b), we 
have 



n 3 (k) = n 2 (k). 



J 2 (k) = V R (k) 



l + 2xA 2 (k) 
A(k) - 29xv R (k)' 
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Inserting this into (5. Id), we obtain 



A 2 (k) = 



2x 



- 1 



A(k) - 29xv R (k) 
[A(k) - 9xv R (k)f 



Now using (5.1a) and the identities mentioned above, 
we find 



A(k) - 29xv R (k) 
[A(k) - 0xv R (k)]* 



= i + ir B (k) + i(r B oA)(k) 



A-2 



A - 28xv R 
[A-9xv R }< 



(k). 



(5.3) 



When 9 = 0, this reduces to (3.11). So far we have not 
made any assumptions on the magnitude of 9. For small 
9, one may consider the lowest order perturbation in the 
structure factor as compared to the pure solution. Note 
that, up to 0(9), the equation takes the sam e form as in 
the 9 = case, except that r#(k) in (5J2) contains an 
additional factor of (1 — 9). This means that for 9 <C 1, 
the solution for given 9 and x is the same as the pure 
solution, denoted by A , at x(l — 9); 



A(k;x,0)ciAo(k;s(l-0)), 9 < 1. 



(5.4) 




FIG. 11. The height Ai of the first peak as a function of 
temperature. The filled circles, triangles, filled squares and 
open squares correspond to the solutions at 8 = 0, 0.1, 0.2 and 
0.3, respectively. The solid lines are linear fits. 



For given x and 9, olt is determined via (2.13). If we 
denote by ) the re lati on between ax and x for 9 = 0, 
then from (2.13) and (5.4) we deduce that to 0(9) 



ar(x) 



a° T (x(l - 9)) - -(2y/x + a° T (x)). 



(5.5) 



We note that, for 9 = 0, a^(x) = (1 — x(3a(x))/\/~x ~ 
— \/x~I3a as ar — » — oo, where the generalized Abrikosov 
ratio changes from 2 at high temperatures to some low 



temperature limit, Pa- In the pure case, we foundi that 
the length scale grows as ~ |oy(x)| ~ s/x~Pa(x)- Using 
(|5.4|) and (5.5), one can derive that 



R c ~ p A J(l=ffjx~ ~ (1 - 9(\ + i))|a T | (5.6) 

2 PA 

as ar — * — oo. As mentioned before, this behavior is 
already captured in our numerical solution (see Fig. fLl]) . 



VI. GENERALIZATION TO THREE 
DIMENSIONS 

In this section, we apply the parquet resummation 
method to a three-dimensional vortex liquid system. We 
shall demonstrate that, while the generalization of the 
parquet equations to three dimensions is straightforward, 
it is virtually impossible to solve the equations using the 
same numerical technique as in the 2D case. We shall 
then discuss a possible analytic approach to the prob- 
lem. 

We consider only the pure case for simplicity. The 
order parameter ^>(x,y,r±) in 3D depends on the coor- 
dinate r± perpendicular to the (x, y) plane. In the lowest 
Landau level approximation, the order parameter takes 
the form 

$(x,y,r±) = <p(z,r±)exp(- — \z\ 2 ). 



The GL free energy analogous to (2.3) is given by 



F[<f>] = J dr±dz*dz((\d±(t>\ 2 + a H \<l>\ 2 )e- 







+ |exp(-^|zn|^,r ± )rj, 

where d± = d/dr±. The renormalized propagator can be 
written as 



G(C,z; q ) = fg(q)exp(^-Cz), 



(6.1) 



where q is the Fourier momentum corresponding to r±. 
For the bare propagator, Qo(q) = (q 2 + o^) -1 . But, 
the renormalization drives Q (q) into a general function. 
This is in contrast to the two-dimensional case where the 
propagator is completely described by one parameter a R . 

Similarly, the renormalized four-point vertex carrys 
four momenta, qi,q2,q3, and q^, in addition to k describ- 
ing the correlation in (x, y) plane. Among the four mo- 
menta only three will be independent due to the momen- 
tum conservation (q\ +q2 = 93 + 94). Using, for example, 
new variables, s = q\ + q^, t = qi — 93 and u = q\ — 94, 
one can describe the renormalized four-point vertex func- 
tion in terms of T(s,t,u;k). The parquet equations for 
r follow from Fig. [|. By evaluating the diagrams in the 
Fig. |B| using the momentum conservation on each vertex, 
one finds that 
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IIi(s,i,u;k) 
U 2 (s,t,u;k) 
n 3 (s,t, u;k) 



->1 



dq 
2^ 



Q(q,s) 



. ( s + t + u s 
Ai(s, 5 g,g 



t — u 



).r( 



s — i + it 
s,9 s A 



2„ I ^-Q(q,t)A 2 (q+- ^ 



g;k r g 



2 

t — u 



s + t — us + t + u \ _ / s 
A 3 9+ s i n *T g + — 



t — u 



s + t — u 



,t,q- 





2 




s 


+ i - 






2 




s 


-i + 




2 



(k) 



(k), 



where r\ = (3^ 2 /2n and 0(g, g') = G{q)G{q — q') and o 
and * operate on the two-dimensional momentum k as 
before. The remaining parquet equations take the same 
form as in the 2D case: 

3 

r(s, t, u; k) = R(s, t, u; k) + ^ IIj(s, i, u; k) 

i=l 

Aj(s,t,u;k) = r(s,t,u;k) - Hi(s,t,u;k), 

where R accounts for the totally irreducible parts, which 
is set to 

R{s, t, u; k) ~ exp(~k 2 /2M 2 ) 

independent of s,t,u in the parquet approximation. 

The first difhculty one faces when one tries to solve 
the above parquet equations as in the 2D case is the fact 
that the propagator is given as an unknown function in 
this case. In the 2D problem, we were able to deter- 
mine the unknown constant an (or x) self-consistently 
using Dyson equation. In this case, we have a func- 
tional self-consistent equation for the self-energy S(g) = 
G~ l {q) - Gq 1 ^) as follows: 



xg(q')g(q")g(q + q'-q")r(q + q',q 



d 2 k dq'dq" e _ k 2 



(27T) 



// // 

q ,q 



12^ 
g';k) 



Even if one uses this self-consistent relation, the par- 
quet equations are coupled integral equations for the 
function T of four independent variables (s, t, u and |k|). 
Storing all the data for T and for all the irreducible 
parts, Aj in the numerical calculation will be a formidable 
task. Furthermore, the integrations in the equations be- 
come multiple sums and therefore the previous numerical 
method of inverting a matrix become more cumbersome. 
One may have to resort to a direct iteration, which con- 
verges slower than the numerical inversion of matrices. 

These numerical difficulties, however, should not dis- 
courage one from attempting to get some nonperturba- 
tive information from the parquet equations. In 2D, 
the parquet equations seem to be a minimal set of 
equations that predicts the growing translational order 
in (x, y) plane, whose length scale is characterized by 
R c . Thus, we expect that the above 3D parquet equa- 
tions contain among other things nonperturbative infor- 
mation on various length scales characterizing the low 
temperature regime of the system. In addition to R c , 
the 3D system may have growing length scales for the 



correlation in the r± direction. We denote by £l the 
length scale arising from the propagator G and by L c the 
one associated with the 3D stru ctu re factor S(k, q) = 
(2ir/fi 2 )ex-p(k 2 /2n 2 )x(k,q) (see (p~4j)), where 



S(k,q) = J ^-G( q ')G(q + q ')-2n J ^G(q')G(q") 
X G(q + q')G(q + q")T{q + q' + q" , 9, q' - 4' \ k). 

If one uses in the parquet equations some kinds of low- 
temperature asymptotic forms for Q{q) and 5(k, q) which 
are characterized by the length scales, L c and R c , 
one might be able to extract nonperturbative relations 
for these length scales and the temperature. We have 
tried a simple ansatz where G and S are represented as 
delta function-like sharp peaks around q = and k = G 
(RLV). The length scales are set to be equal to the in- 
verse width of the corresponding peaks. But, we found 
that these forms are too simplified to produce sufficient 
amount of information on the temperture dependences 
of each length scale. We believe that more sophisticated 
asymptotic forms arc needed for the analysis of the 3D 
parquet equations. This point will be discussed again in 
the next section. 



VII. DISCUSSION AND SUMMARY 

In this paper, we applied the parquet resumma- 
tion method to the 2D vortex liquid with and with- 
out quenched impurities and to the 3D vortex liquid 
in the absence of disorder. In the 2D system, we were 
able to solve the parquet equations numerically and find 
the length scale R c . The temperature dependence of 
R c was also obtained. In the pure 2D vortex liquid, 
we found in the previous paperEl that the asymptotic 
forms, A(k) ~ 27r^ 2 J2g^o ( ^ 2 * > (k ~ G), or equivalently 
f R (k) - - 7r^ 2 X; G ^ (2) (k - G)] solves the par- 

quet equations if the inverse width of the delta-function 
peaks, or the length scale R c behaves like s/x ~ \ar\ as 
c\T — * —00. In the disordered case, one can obtain similar 
low-temperature asymptotic forms for /i?(k) and w_R(k) 
in terms of delta functions around k ~ G. One can then 
obtain expressions for R c similar to (p 7 ^). These simple 
asymptotic forms, however, have some limitations. First 
of all, the amplitude of each delta-function peak cannot 
be determined from the parquet equations. It is also not 
clear whether one can use these representations in the 
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strong disorder regime. In the previous section, we ar- 
gued that an ansatz like S(k, q) ~ 5(q) J2g^o & s - — G) 
was not enough to produce any useful nonperturbative 
information for the three dimensional solution. We be- 
lieve that finding appropriate low temperature asymp- 
totic forms for the vertex functions is a first step towards 
a more fruitful use of t he par quet equations. To this end, 
a simple equation like ( 3.11 ) might be useful. 

The effect of random impurities on the mixed state of 
a type-II superconductor is usually described by the so- 
called Larkin-Ovchinnikov arguments, which states that 
for a dimension d < 4 the long-range crystalline order of 
the mixed state is destroyed by a weak disorder. For thin 
films (d = 2), one can estimate the Larkin length scale R c 
over which a short-range order persists ascl R c ~ |ar|/A. 
Although the |ar (-dependence of R c is the same as the 
present result, we note that there is a basic difference 
between our result and the LO type argument. Since the 
LO argument starts from a perfect crystalline state where 
A = 0, R c = oo, a very small amount of disorder makes 
an abrupt change as can be seen from 1 / VA- dependence 
of R c . Within the parquet approximation, there is no 2D 
crystalline phase and the pure system has a finite R c at 
any finite temperature. As we have seen in the previous 
sections, the effect of small disorder is represented as a 
smooth decrease of this length scale. 



In this paper, we have concentrated on a rotationally 
symmetric liquid state where the vertex functions depend 
only on the magnitude of k. But, for example, in the 
presence of magnetoelastic interactions between vortices, 
one must consider general k-dependent vertex functions 
in the parquet equations. We recall that the input pa- 
rameters of the parquet equations are the temperature 
and the bare vertex functions, fs and vb- It would be 
interesting to study a situation where the bare vertex 
functions possess a nontrivial k-dependence. 

One could also generalize the present parquet resum- 
mation technique to study the dynamics of type-II su- 
perconductors. We expect that the parquet equations 
for the time-dependent Ginzburg Landau theory can be 
obtained without much difficulty. But, the vertex func- 
tions will depend on the frequency as well as the usual 
momenta. Thus, one faces similar numerical difficulty to 
the 3D case. la addition, one has to be careful in taking 
the LLL limitca. 
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APPENDIX: EVALUATION OF DIAGRAMS 



Here we give some examples of evaluating the Feynman diagrams involved in the present work. Similar analysis 
can be found in the Appendix of Ref. |2l]. We consider only the pure case (A = 0) for simplicity. The results can easily 
be generalized to the disordered case carrying the replica indices. Feynman diagrams are constructed by connecting 
the propagator lines represented by (2.6) with the quartic vertices, which can be obtained from (2.7) as (see Fig. |T2j) 







dkdk* 



r(k)e 



exp 



(V(c-c') + (c*-c'*)&) 



By joining four propagator lines (starting from z* and and ending at z% and 24) with this vertex, we have 

2 



-4(f)(f) 4 4 
z 27T off 



dkdk* 



-r(k)c 



fe| J /V 



%-^(ICl| 2 + IC2| 2 )/2 



x exp 



(*rci + 4b + <r*3 + C2Z4) - \ (V(d - co + (cr - c 2 *)fc) 



After performing the Gaussian integrals using 



P m 

J HdQdQeM-AQMnQ - a*Q - = (^) m (det M)" 1 exp [/x 2 (a* Mr ^ 



(Al) 



one obtains (2.9). 

Let us evaluate the one-loop diagrams and derive ( |3.6| ). We consider only the first diagram of Fig. || and thus derive 
( |3.6a ). The remaining two diagrams can be evaluated in the same manner. The first one- loop diagram gives 



2| 2 )/2M 2 



JJ dtjdCj exp 
3=1 



2 4 

3=1 
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x exp 
x exp 



(z*Ci + Z2C2 + (3Z3 + CU4 + C1C3 + C2C4) 

\ (k(Ci - c 2 ) + (cr - c 2 *)fci + k (c 3 - a) + (c 3 * - en**) 



Integrating over and Q using ( |A1[ ) and changing the variables + k 2 — > k, and k x — > p, we obtain 



( M ) 2 ( M ) 2 c M 2 (^r^3+^^4)/2 

2 27T 



^A l( p)r(k-p )e -^^ 



exp 



(fc*(z 3 - Z 4 ) + (z* - Z*)fc) + (fc> - p*fc) 



Comparing this with the general expression ( |2.S| ), we find that the one- loop contribution from this diagram to the 
vertex function is 

/3/x 2 f dpdp 



2na R 



J ^ A i(p) r (k-p)exp(^(rp-p*fc)), 



which is equal to the right hand side of (3.6a). 



-f r N 

FIG. 12. Diagrammatic representation of the propagator and the four-point vertex. 
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